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We discuss the problem of competition between a superconducting (SC) ordered state with a 
charge density wave (CDW) state in stripe phases of high T^ superconductors. We consider an 
effective model for each stripe motivated by studies of spin-gapped electronic ladder systems. We 
analyze the problem of dimensional crossover arising from inter-stripe SC and CDW couplings using 
non-Abelian bosonization and renormalization group (RG) arguments to derive an effective 0(4)- 
symmetric nonlinear cr-model in D = 2 -I- 1 for the case of when both inter-stripe couplings are 
of equal magnitude as well as equally RG relevant. By studying the effects of various symmetry 
lowering perturbations, we determine the structure of the phase diagram and show that, in general, 
it has a broad regime in which both orders coexist. The quantum and thermal critical behavior is 
discussed in detail, and the phase coexistence region is found to end at associated T = as well 
as T > tetracritical points. The possible role of hedgehog topological excitations of the theory 
is considered and argued to be RG irrelevant at the spatially anisotropic higher dimensional low- 
energy fixed point theory. Our results are also relevant to the case of competing Neel and valence 
bond solid (VBS) orders in quantum magnets on 2D isotropic square as well as rectangular lattices 
interacting via nearest neighbor Heisenberg exchange interactions. 

PACS numbers: 74.81.-g,74.40.Kb,74.62.-c 



I. INTRODUCTION 

Over the past few years substantial evidence, both 
experimental and theoretical, has accumulated showing 
that the pseudogap regime of high Tc superconductors 
(HTSC) is related with to existence of charge and/or 
spin ordered states, stripe phases that break translation 
symmetry in the conducting layers, as well as uniform 
anisotropic nematic states (for reviews see seei"— and ref- 
erences therein). In addition, homogeneous phases with 
"hidden order" involving some degree of time reversal 
symmetry breaking have also been proposed for the pseu- 
dogap phase*^'^. An important question in this field is to 
determine if these ordered phases either compete with or 
are part of the mechanism of high Tc . 

In principle, superconducting (SC) and charge/spin 
stripe orders can coexist or compete. For instance, in 
a well developed SC (gapped) coupled to a weak charge 
density wave (CDW), the CDW has little effect. Even 
deep in a d-wave SC, weak stripe order (a unidirectional 
CDW) does not generally affect the low energy physics 
of the SC quasiparticles as its wave vector generally does 
not span the nodal points^'^. Yet, in a weak coupling 
BCS approach to both orders, they usually (but not al- 
ways) compete with each other, and the enhancement 
of one type of order normally suppresses the other, and 
one has "competing orders" (see e.g. Ref.[9] and refer- 
ences therein). However, recent work on stripe phases 
has provided a different perspective on this question by 
showing that in a strongly correlated system inhomoge- 
neous phases not only are unavoidablel*^ but may also be 
part of the mechanism of SC and, moreover, that there 
exists an optimal degree of inhomogeneity at which the 



Tc is largestiii^. 

In this work we consider the problem of the interplay 
between SC and stripe phases of HTSC and formulate an 
effective field theory of these competing/ coexisting or- 
ders. We will use as a starting point a strong coupling 
picture of a stripe phase in which this 2D ordered state 
is represented as a quasi-lD system, an array of doped 
Hubbard-type ladders each in a Luther-Emery liquid, a 
phase with a finite, and typically large, spin gap (de- 
noted by Ag) ^^d^ . There is extensive numerical evidence 
for this assumption based, primarily, on density matrix 
renormalization group (DMRG) calculations of Hubbard 
and t — J models on up to seven leg laddersi^"— . There 
is also supporting analytic evidence for spin gap phases 
in doped ladder s^^d Z, although the relation between the 
parameters of the effective (bosonized) field theory thus 
derived and microscopic models can only be reliably de- 
termined from numerics. In the low-energy regime these 
systems exhibit spin-charge separation. Furthermore, 
for a broad range of doping, < a; < 0.3, ladders 
with only repulsive interactions have a spin gap but do 
not have a charge gap (provided a; > 0) (for a review 
see^.) The states on the ladders represent the local 
high energy physics from which the low energy 2D states 
arise. (These models can also be used to describe ladder 
materialsi^) . Thus, the resulting SC state is not the re- 
sult of pairing in a system with preexisting quasiparticles. 
Instead, the quasiparticles of the SC state are emergent 
low energy excitations. 

Based on this strong coupling picture of the stripe 
phase we use bosonization methods to determine the 
structure (and symmetries) of the effective low energy 
physics at scales well below the spin gap and construct 



an effective field theory for the resulting 2D state, a 
problem that was has not been addressed before. The 
power of this approach resides on the fact that allows 
the treatment of the competing effects of the tendency 
of the system become a 2D superconductor and a 2D 
charge ordered phase (a CDW). It turns out that, for 
a range of parameters, the system naturally has an ap- 
proximate (enlarged) symmetry, that includes both the 
SC and stripe orders, that has far reaching consequences 
on the structure of the phase diagram. In particular, 
this approximate symmetry allows us to establish the 
existence of an intermediate phase in the 2D system in 
which CDW and SC orders coexist rather than being sep- 
arated by a first order transition. This general question 
was discussed earlier on using inter-ladder renormaliza- 
tion group (RG) methods^^''^*' and inter-chain mean field 
theory (ICMFT)^, suggesting (in both cases) that the 
transition (at T = 0) is first order and that there is an 
associated bicritical point. In this paper we will show 
that, instead, there is a phase in which SC and CDW 
coexist and that there is a tetracritical point (both at 
T = and at r > 0.) 

This work is is organized as follows. In Section |TT] 
we present a summary of the stripe model and use it 
to determine the structure of the effective field theory. 
The effective field theory is derived in Section Hill The 
quantum and and thermal (classical) critical behaviors 
are discussed in Section ITVl We close with a summary of 
results and extensions in Section [V] 



II. ANALYSIS OF THE STRIPE MODEL 

We consider an array of weakly coupled spin-gapped 
ID systems with a gapless charge sector ("ICSO")^. We 
will follow the construction of Refs. [Il|[l9| . The low 
energy Hamiltonian for an array of decoupled ladders has 
the bosonized form 



leading relevant inter-stripe operators are the Josephson 
and CDW interactionsii-i^iiSi^a 
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(2.1) 

where i is the stripe label, Kc is the charge Luttinger 
parameter of each ladder; (pi is the phase field of the 
CDW along each ladder and 6i are the SC phase fields 
of each ladder. We will ignore other gapped degrees of 
freedom of the ladder as we are only interested in the 
physics below the spin gap Ag. 

The fields (pi and 6i are dual to each other, and satisfy 
the equal-time commutation relations [(j)i{x),dxOj{x')] = 



idijS{x 



The SC and CDW order parameters on 



each ladder in a spin gap phase are given by Osci^) ~ 
exp{i\/2n6{x)) and OcDwix) ^ exp(iV27r0(x)), respec- 
tively. 

DMRG studies of hole-doped ladder systems^ 
showed that K^ is a function of doping x with Kc — >■ 2 as 
X — >■ to Kf. — >■ 1/2 a.s X —i' Xc — 0.3. In the presence of 
a spin gap the single particle hopping is irrelevant. The 
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(2.2) 

Notice that we have used the seemingly wrong sign for 
the Josephson coupling (as it induces a tt phase shift be- 
tween neighboring ladders). However this is the natural 
sign for the coupling between the CDW order param- 
eters on neighboring ladders since they want to couple 
with a TT phase shift so as to lower the effects of the 
repulsive interaction. We will refer to this as an "anti- 
ferromagnetic" sign. However, since the two phase fields 
(j)j and 6j are conjugate to each other (as required by 
their commutation relations), it is possible to carry out 
a unitary transformation on the SC phase fields 6j on 
every other ladder (induced by the action of unitary op- 
erators involving only the CDW phase fields of the same 
ladder) and map the problem to one in which both in- 
teractions have a "ferromagnetic" sign, and both phases 
want to lock-in without a tt phase shift. Thus, in this 
simple state in which there is a spin gap on every ladder, 
the state associated with the sign alternation of a phase- 
density-wav o^^'^^ is equivalent to one without the phase 
alternation. In contrast if both interactions had a "ferro- 
magnetic" sign from the outset (i.e. if the phases wanted 
to lock without a tt phase shift) it is no longer possible to 
map the system to one in which both are "antiferromag- 
netic". We will see below that this makes an important 
difference to the form of the effective field theory. 

In the spin gap regime, Q < x < xc — 0.3, the scal- 
ing dimensions of the CDW and Josephson couplings of 
Eq. (|2.2p is A^c = 1/^c and l^cDW = Kc respectively. 
Hence, both couplings are relevant since their scaling di- 
mensions satisfy l^sc,CDW < 2, with SC being more 
relevant than CDW for x < Q.l (where Kc — 1) and 
conversely for 0.1 < a; < 0.3 (where As -^ 0). Then, 
a perturbative RG analysis of these two inter-stripe in- 
teractions, valid for Jsc and Jcdw sufficiently small, 
shows that for < a: < 0.1 the system fiows to a 2D 
(anisotropic) SC (with subdominant CDW correlations) 
whereas for 0.1 % x "^ 0.3 it flows to an (anisotropic) 
crystal (a bidirectional CDW, with subdominant SC cor- 
relations) suggesting a direct first order transition at 
X ~ 0.1^^'^^. We will only consider nearest- neighbor 
inter-stripe interactions, as couplings for further away 
stripes become small quickly (although they are just as 
relevant.) 

Inter-chain mean field theory (ICMFT)i^i2£i2Z, also 
found a direct SC-CDW first order transitional. ICMFT 
yields a good qualitative picture of the dimensional 
crossover and, indeed, it has the correct asymp totic scal- 
ing as Jsc/CDW -^ (see, e.g. Refs.[ll|[illlil), pre- 
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TcDW « JcDW^"^^ respectively. Within ICMFT, close 
enough to the critical temperatures Tsc and Tcdw, 
the order parameter is small and follows the critical 

behaviors2i.29 ^.^(T) oc Asc(O) (l - ^) ''^, and 



^CDw{T) oc AcDwlO) ( 



1 - 



l/(4-ifc) 



These re- 



suits are consistent with those of Ref . (13| for the case of 
Kc = 1/2 (weakly coupled Luther- Emery stripes), which 
predicts a 2D SC phase with BCS- type scaling. The 
problem of the competition of SC and CDW ordering was 
examined in considerable detail in Ref. '21] at the level of 
ICMFT, supplemented by a random phase approxima- 
tion (RPA) analysis of fluctuation effects, and concluded 
that, at this level of approximation, there is a direct first 
order transition from SC to CDW order at T = . 

However, the problem of whether in the 2D system 
there is a direct SC-CDW first order transition (and an 
associated bicritical point) or a phase in which both order 
parameters coexist (controlled by a tetracritical point) 
cannot be addressed at the level of ICMFT or within 
the quasi- ID regime, and requires an RG analysis of the 
2D problem. In answering this question, we take an al- 
ternate route, and derive an effective field theory of the 
2D system, and study it as a problem in (quantum and 
thermal) critical phenomena. 

In the regime where the SC/CDW competition is 
strongest, i.e. near the self-dual point (SDP) at K^ = 1, 
both the SC and CDW interactions have the same scal- 
ing dimension Age = ^cdw = 1, and hence must be 
treated on an equal footing. Under duality, 9i o 4>i 
and Kc O K~^ , and at Kc = 1 the Hamiltonian 
H = i/intra + -fainter of the 2D systcm is exactly self- 
dual (if the inter-stripe couplings are equal). Moreover, 
at Kc = 1 the system has a (dynamical) SU{2) symme- 
try, a feature reproduced by ICMFT—. The chiral fields 
on each ladder 0^;/^ = (6* ± (/')/2 can be used to define 
three chiral current operators 
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J± ^ ~ exp (±i2V2^,pR^L) (2.3) 



each with scaling dimension 1, which generate an SU {2)i 
Kac-Moody algebra independently for the right and left 
moving degrees of freedom^. We will see below that 
at Kc — 1 SU{2)i^ X SU{2)l is also a symmetry of the 
full 2D Hamiltonian, Eqs. ((2l1) and ((2?^ . provided the 
inter-stripe couplings are equal. This symmetry is bro- 
ken down to U{l)!i X C/(1)l both away from Kc — 1 
and by unequal couplings. The total interacting Hamil- 
tonian, H ~ -ffintra + -Winter is invariant under global 
U{1)b. X U{^)l (chiral) symmetries. In the absence of 
inter-stripe couplings these global symmetries turn into 
sliding symmetriesi^iSS. We note that the Hamiltonian 
in Eq. (|2.1I) at Kc = 1 also represents the continuum limit 
of quantum Heisenberg antiferromagnetic chains. With 
some caveats discussed below, our results will also be rel- 
evant to the case of an array of such spin chains weakly 



III. EFFECTIVE FIELD THEORY 

Let us now derive first an effective field theory 
in 2 -f 1 dimensions in the regime in which these SC 
and CDW maximally compete with each other. Hence 
the effective field theory will have, to zeroth-order, an 
SU{2) X SU{2) symmetry that will be explicitly broken 
down to U{1) X U{1). In spirit, this approach is reminis- 
cent to the phenomenological 5'0(5) theory^ although 
the microscopic origin of the symmetry is quite differ- 
ent. Provided Kc ~ 1, it is possible to treat the symme- 
try breaking terms in a controllable fashion. To do this 
we begin by using the fact that, at Kc = 1, the effec- 
tive Hamiltonian of each ladder (j = 1, . . . , iV), Eq. (|2.1|) . 
has a well known equivalent representation in terms of 
an SU {2)i Wess-Zumino-Witten model (WZW) whose 
Hamiltonian density i o^o, 32^33 
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E (■■ jjA^yiBi^) ■■ + ■ jjA^yui^) 



(3.1) 
where the (normal ordered) products involve the right 
and left moving currents Jf^ji{x) and Jj^l{x) (a ~ 1, 2, 3), 
the three generators of two 5*^7(2)1 chiral Kac-Moody al- 
gebras. Away from Kc = 1 there are symmetry-breaking 
{SU{2)ii X SU{2)l -^ U{1)r X U{1)l) intra-stripe cou- 
plings of the form 
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SB 



(3.2) 
with w — > is Kc — ?► 1, as well as inter-stripe (gen- 
erally anisotropic) current couplingslS. As shown in 
Rcfs.[19,20], sufficiently strong and long-ranged inter- 
stripe current-current interactions can turn the inter- 
stripe SC and CDW interactions RG irrelevant and lead 
to the low-energy fixed point of the sliding TLL (or smec- 
tic liquid). We will focus here, however, with the oppo- 
site case, i.e., that of dominant SC and CDW inter-stripe 
interactions. 

The WZW model is a non-linear sigma model 
(NLSM) in 1 -I- 1 dimensions whose degree of freedom is 
a field g{x) that takes values on a (compact) Lie group, 
SU{2) in this case. The action of the WZW models for 
each stripe is^ 
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5wzw[5] = Y^ I d^xiT{d^g^d,,g) 
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(3.3) 

where k is the level of the Kac-Moody algebra; here we 
are interested in the case SU {2)i and A; = 1. The second 
term in Eq. (|3.3l) is the WZW term, where B denotes a 



3D solid sphere whose boundary 5^ is 1 + 1 dimensional 
space-time. The S'C/(2)-field g and (</>, 0) fields are related 

by 
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The (combined) inter-stripe terms in the Hamiltonian 
density now take the form 

"Hintcr == -J+^tr(5](a;)5j+i(x)) 

- J-^'^A9]^^y^93^i{x)^z) (3.5) 

3 

where Uz is the diagonal Pauli matrix, and we have de- 
fined the couphngs J± = Jsc ± Jcdw- Notice that 
for J71 = the 2 + 1-dimensional system still enjoys 
an SU{2) x SU{2) symmetry, which is broken down 
to [/(I) X [/(I) if J- ^ 0. In addition, away from 
Kc — 1, the intra-stripe Hamiltonian has additional 
terms of the form tr{dfj,g^azd^gaz) that also break the 
symmetry down to t/(l) x U{1). 

Depending on the signs of J7± , the inter-stripe cou- 
plings favor ordered phases in 2 -I- 1 dimensions in which 
the SU{2)-vahied matrix field g is either uniform across 
the system or changes sign (i.e be staggered by an ele- 
ment of the Z2 center of the group SU{2)) from one stripe 
to the next. In a (relatively) recent paper, Senthil and 
Fisher'^^ discussed the behavior of an array of antifer- 
romagnetic quantum Heisenberg model in the quasi-lD 
regime and proposed a description of that system which 
is similar in spirit to the one we use here for the stripe 
state in the spin gap regime. One important difference is 
that in the Senthil-Fisher the SU{2) group is associated 
with the spins which, in the ID limit, have quasi-long 
range antifcrromagnetic commensurate order with wave 
vector TT. As a result, in Ref. [34!] the interchain coupling 
has the form ii{gigij^i) + h.c. that breaks the symmetry 
to a single global SU{2) symmetry. This coupling favors 
a staggered order in which g alternates with g^ between 
neighboring chains. In the problem we discuss here this 
form of coupling is not allowed by gauge invariance (in 
the case of the SC order) and by translation invariance 
(in the case of CDW order). As we will see below this 
leads to some differences in the structure of the effective 
field theory. 

The derivation of the effective field theory in 2 + 1 
dimensions proceeds in two stages, e.g. along the lines 
discussed early on by Affleck and Halperin^^. In the 
quasi- ID limit the leading relevant operators are the 
inter-stripe couplings shown in Eq. p.Sp . In the ID limit 
at the SU{2) invariant system (with Kc — 1), in the 
absence of the inter-stripe interactions, the low energy 
physics is controlled by the WZW fixed poinli^, an in- 
frared stable fixed point at a finite value of the NLSM 
coupling constant. At the WZW fixed point the in- 
ter stripe operators have scaling dimension 1, strongly 



destabilize the fixed point of the WZW model, and fa- 
vor 2 -(- 1-dimensional ordered states at which the field 
g picks up a non-vanishing expectation value. In the 
presence of the "internal anisotropy" terms that break 
SU{2) X SU{2) ^ Uil) X U{1) the resuhing ordered 
phases correspond to 2D SC, 2D CDW and possibly co- 
existence phases (which is the main question we address 
here). 

After this initial step of quasi- ID renormalization, 
the system becomes coarse-grained and flows to a full 
2 + 1-dimensional theory with the proper symmetries, 
with an effective field theory of the form of a relativistic- 
like 2 -I- 1-dimensional NLSM. The effective theory is gen- 
erally spatially anisotropic and also has terms "internal" 
anisotropy terms. However, the infrared unstable (quan- 
tum critical) fixed point of the 2 + 1-dimensional NLSM 
occurs at a finite value of the NLSM effective coupling 
constant and, hence, cannot be accessed perturbatively 
from the quasi-lD regime (at least not in a controllable 
way). Thus, derivations of the effective theory based on 
naive gradient expansions of the quasi- ID Hamiltonian 
must be regarded as being only suggestive (at best) of 
the structure of the effective field theory near the quan- 
tum critical point. Nevertheless it is possible to use the 
powerful constraints of locality and symmetry to write 
down the structure of the effective field theory. This is 
the approach we use here. With one significant caveat 
concerning the role of topology (that we will discuss be- 
low), the spatial anisotropics discussed above lead to re- 
dundant operators whose effects can be taken into ac- 
count by a suitable rescaling of the spatial coordinates 
and time. In contrast, the internal anisotropy terms play 
a key role. In what follows we will work in imaginary 
time. 

For reasons of clarity we find it useful to represent the 
5t/(2)-valued matrix field g in terms of a four-component 
vector iV" (a = 0,1,2,3) that takes values on the three- 
sphere S^, i.e. A^° = ^tT{ga°'), where we have used the 
basis of Hermitian 2x2 matrices, a^ — I and the three 
Pauli matrices for a = 1,2,3. The four-component field 
N'^ satisfies the constraint N"^ = 1, and as such takes 
values on S^. 

Given these symmetry requirements, the only al- 
lowed effective action of the effective field theory in 2 -I- 1 
dimensions is that of an 0(4) NLSM. Ignoring for the mo- 
ment spatial anisotropics and internal anisotropy terms, 
the effective Lagrangian is 
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'2g^ 
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(3.6) 



where fi — 0, 1,2 label time and the two spatial coordi- 
nates respectively. The bare value of the effective cou- 
pling constant go of the NLSM in 2 -I- 1 dimensions has 
units of (length)^^. Its value is essentially given by the 
geometric mean of the (suitably dimensionalized) cou- 
pling constant of the WZW model, l/(47rao) (where ag 
is the stripe spacing) and the inter-stripe coupling J7+. 
This value of go will be substantially modified by renor- 



malization effects. 

The competition/coexistence of SC and CDW orders 
can be now be discussed by considering the effects of 
operators that break the large SU{2) x SU{2) symmetry 
down to U{l)xU{l). The symmetry lowering arises from 
(i) a finite asymmetry in the bare values of the SC and 
CDW couplings, Jsc/cdw 7^ 0, and (ii) a departure 
from the SDP as Kc — 1 7^ 0. Up to spatially anisotropic 
gradient terms that we will omit for now, we obtain the 
following effective Lagrangian 

C,s[N] = -^{d^Nf - W d^NaOabd^Nb - h NaOabNb 

(3.7) 

where a, & = 0, 1, 2, 3. In Ea. (l3.7p the second term with 
coupling constant w ^ K^ — l parametrizes the departure 
from the SDP, the third term with coupling constant h ^ 
J71 describes the unequal couplings between CDW and 
SC order parameters, and O is a fixed matrix that breaks 
SU{2) X SU{2) -^ f/(l) X C/(l). 

There is still one more operator that can be part 
of the effective action that needs to be considered: the 
topological charge Qtop defined by 
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(3.8) 

Qtop[-^] is an integer, a topological invariant that classi- 
fies the non-trivial maps of the 2 + 1-dimensional Eu- 
clidean space-time (compactified to S^) to the target 
manifold of the 0(4) NLSM which is also isomorphic to 
S^. In the language of homotopy theory these maps are 
classified by n3(S''^) ~ Z. Since the system at hand is 
time-reversal invariant, the only allowed topological term 
in the effective (Euclidean) action must have the form 
'S'top = ipT^Q top[N], where p G Z. Since Qtop[N] is also 
an integer, such a topological term has an effect only for 
p odd. 



IV. CLASSICAL AND QUANTUM CRITICAL 
BEHAVIOR 



A. r = o 

We begin the discussion of quantum criticality by 
considering the role of the topological term. Terms of this 
type play a crucial role in ID spin-1/2 quantum Heisen- 
berg antiferromagnets'^^, also described by an SU{2)i 
WZW model^. Extensions to higher dimensional an- 
tiferromagnets have also been proposed before and ex- 
amined in some detail^>^"— , but never in the context of 
stripe phases in superconducting systems. Particularly 
relevant here is the work of Senthil and Fisher'^^ who 
discussed this same effective field theory in the context 
of the quantum phase transition from a 2D Neel antifer- 
romagnetic state to a four-fold degenerate valence bond 



solid (VBS) state on a square lattice, whic has been con- 
jectured to be controlled by a deconfined quantum critical 
(DQC) point^. The topological excitations that carry a 
non-trivial winding number Q[N] are monopole ("hedge- 
hog") configurations in 3D Euclidean space- time, the 
instantons of this theory. Contrary to their 2D cousins, 
the Euclidean action of instantons of 3D NLSMs is lin- 
early divergent and hence are suppressed throughout the 
ordered phase. Nevertheless, within the DQC scenario 
they are still argued to play a key role both at the quan- 
tum critical point and in the quantum disordered phase 
(which becomes a topological phase). However, the anal- 
ysis of Ref.jSJ] shows that the 3D DQC fixed point is 
unstable to the effects of spatially anisotropic perturba- 
tions, such as the ones we have in this theory, and be- 
comes inaccessible. The resulting effective field theory of 
our system is then the anisotropic NLSM of Eq. (|3.7p but 
without the topological term which becomes an irrelevant 
operator at the accessible fixed point. 

Having determined the form of the effective field the- 
ory we can proceed to study its quantum and classical 
critical behaviors. The analysis of the phase diagram 
can now be determined using well established methods 
of (classical and quantum) critical phenomena and, in 
particular, of multicritical phenomena. At this new, 
more conventional, fixed point spatial anisotropics be- 
come mild redundant operators and their presence no 
longer affect the critical behavior. We will see now that 
the effective action of Eq. p.7p will allow us to find a 
answer to the problem of SC/CDW competition vs. co- 
existence, both at r = and at T > 0. A summary of 
the results is presented in the phase diagrams shown in 
FiglDandFigia 

We thus have a 2 -t- 1-dimensional NLSM with sym- 
metry breaking fields. By power counting we see that the 
third term in Eg. p. 71) is the most relevant perturbation. 
From this point of view the problem of the CDW/SC 
competition is conceptually similar to other problem in 
which there is a partial breaking of the internal symme- 
try, e.g. the spin flop transition in magnets. This point 
of view will give us the solution of the problem. The crit- 
ical behavior obeyed by this system must be approached 
either by means of a) a 2-|-e expansion of the NLSM (here 
2 means 1 -f 1 Euclidean space-time dimensions), or b) 
a 4 — e expansion of the associated Landau-Ginzburg- 
Wilson (LGW) theory, usually known as cj)'^ field theory 
(again, here 4 means 3 + 1 Euclidean space-time dimen- 
sions) . Although we already have the problem expressed 
as a NLSM we will have to use the 4 — e expansion ap- 
proach. The reason is that the 4 — e expansion of the cj)'^ 
has a property known as Borel-summability which allows 
for an accurate determination of its critical exponents 
directly in D = 3 Euclidean space-time dimensions (see, 
e.g. Ref . [4l| ) . In contrast, the conceptually simpler 2-f e 
approach has poor convergence properties and has never 
been successfully used to compute exponents in Z? = 3 
even for the simplest NLSM. 

For this reason we will replace the NLSM effective 



action for the field N with 0(4) symmetry by a theory 
with the same symmetries but in which the constraint 
of the NLSM is replaced by a suitable potential. Let us 
denote by N^j, the upper two (CDW) components and Ne 
the lower two (SC) components of the four-component 
scalar field N, respectively. In this form we are describing 
the breaking of 0(4) — ^ 0(2) x 0(2). This is a particular 
case of the breaking of 0{n) down to 0{ni) x 0(712), 
with n ^ rii + n2, that has been studied in detail in the 
literature. It was originally studied to one-loop order in 
the 4 — e expansion by Kosterlitz et al^, and was more 
recently reexamined by Calabrese et al^ who used a five 
loop 4 — e expansion with Borel resummation and were 
able to determine the critical behavior in _D = 3 with 
high precision. 

The resulting LGW Lagrangian (or "free energy den- 
sity" ) has the form 
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2 

4! 



C[N^,Ne] = -id,NA +-(dM +^Nl + ^N; 



r0 



re 



2w 



+^(»i) +^(».r+i»j". 



(4,1) 

As usual, r0 and rg measure the departure from the (clas- 
sical or mean field) critical point; r — [rg + t^)/2 qual- 
itatively plays the role of the coupling constant qq of 
the NLSM, h — {rg — r^)/2 of the symmetry breaking 
field, and ug, u^ and w are four coupling constants that 
parametrize the potential. For u = ug = u^ = w the 
quartic terms have 0(n) symmetry. Notice also that in 
this form the spatial anisotropics can always be absorbed 
by a suitable rescalings of the coordinates and fields. 

Below 4 dimensions, the free field (Gaussian) fixed 
point is always unstable. Kosterlitz et al^ showed that 
theories of these type can describe either a bicritical point 
(the endpoint of a line of direct first order transitions 
between the two phases with order parameters N^ and 
Ng) or a tetracritical point (the endpoint of a region of 
phase coexistence of the two order parameters). They 
also showed that when the bicritical scenario holds, the 
fixed point associated with the critical endpoint has max- 
imal symmetry, 0{n). In the tetracritical scenario, the 
0(n) fixed point is unstable and two possibilities arise 
for the endpoint: a) either it is a decoupled fixed point 
(DFP) at which the 0(ni) x 0(712) theory is effectively 
decoupled (w -> 0), or b) its is a biconical fixed point 
(BFP) at which the 0(ni) x 0(n2) has a non-trivial cou- 
pling w* . 

The one-loop analysis of Kosterlitz et al^^ predicts a 
tetracritical point with a DFP for some sufficiently large 
value of n > 4, and a bicritical point for the spin-flop 
transition 0(3) — > 0(2) x Z2. However, the one-loop 
results are unable to resolve the case of interest here, 
0(4) -^ 0(2) X 0(2). On the other hand, the five-loop 
results of Calabrese et al^ show without ambiguity that 
in D = 3 the 0(4) -^ 0(2) x 0(2) theory has a tetracrit- 
ical point and not a bicritical point as the 0(4) invariant 
fixed point is unstable. However, the five-loop results 



do not have sufficient accuracy to distinguish between a 
tetracritical point controlled by a DFP or by a BFP. Nev- 
ertheless, regardless of this technical issue, this analysis 
predicts the existence of a phase in which the SC and 
CDW order parameters coexist at T = 0. This is shown 
in the schematic phase diagram of FigH] as the phase 
labeled by SC-hCDW. The relative strength of the SC 
and CDW order parameters varies continuously across 
this phase as the magnitude (and sign) of the coupling 
h to the symmetry breaking field varies. (The broken 
line in Fig IT] is not a phase transition and denotes the 
manifold with higher symmetry). All the phase bound- 
aries are in the universality class of the 3D classical XY 
model (i.e. the 0(2) Wilson-Fisher fixed point in _D = 3). 
Whether the tetracritical point is governed by a different 
coupled fixed point (BFP) or a decoupled one (DFP) is 
not presently established. 

We end this subsection by discussing briefly the ex- 
tension of this analysis to a stripe phase in 3D. For sim- 
plicity we will assume that the stripe phase consists of a 
stack of 2D stripe phases and, hence, that the order pa- 
rameter theory is the same as in 2D. Much of the analysis 
above here follows here too. The main difference is that 
the effective field theory is now in D = 4 Euclidean space- 
time dimensions. In this case the 0(4) NLSM does not 
have a topological term to begin with. In D = 4 Eu- 
clidean space-time dimensions the free-field (Gaussian) 
fixed point is marginally stable and the classical (mean 
field) results are correct up to logarithmic corrections to 
scaling. In this case there is clearly a phase in which SC 
and CDW orders coexist. The resulting phase diagram 
in 3D at T = has the same topology as in FigH] 



f/(l) 
Tetra-critical 




FIG. 1: (color online) Schematic phase diagram for a 2D 
stripe phase at T = 0: go is the coupling constant of the 
0(4) NLSM and h is the anisotropy. Here we assumed that 
the tetracritical point is a decoupled fixed point. (See text 
for details.) 




3L>XY 



3D XY model 

(a) T > 0, 30 > 3c 




h 




SC + CDW 

(c) T > 0, 30 < 3c 



FIG. 2: (color online) Schematic T > phase diagrams for a 
2D stripe state as a function of the anisotropy h of the 0(4) 
NLSM for values of the coupling constant go corresponding to 
the lines labelled by a, b and c in Fig[T] See text for details. 



B. r> 

We now turn to the phase diagram at finite temper- 
ature, T > 0. It is not possible to replace the NLSM 
by the LGW theory in 2D at finite temperature, since in 
2D there is a drastic difference in the thermal critical be- 
havior of the system if it is fine-tuned to have the larger 
0(4) symmetry from where it has a lower 0(2) x 0{2) 
symmetry. (Recall that topological terms do not play a 
role in the thermal critical behavior as they always con- 
tain time derivatives of the fields.) (This was already 
emphasized by Carr and TsvelikSi.) In both cases there 
is no long range order in 2D as required by the Mermin- 




4i:»XY 



FIG. 3: (color online) Schematic T > phase diagrams for 
a 3D stripe state as a function of the anisotropy h of the 
0(4) NLSM for values of the coupling constant go < gc- The 
thermal transitions are 3D XY and the T = transitions 
are 4D XY. TCP labels the tetracritical point. See text for 
details. 



Wagner theorem. However, the critical temperature of 
the 0(4)-invariant classical 2D NLSM is zero and it is in 
a classically disordered phase at all T > 0, as shown in 
the schematic phase diagrams of Figl2](a-c). 

Away from the 0(4) symmetric theory, i.e. for h ^ 
0, there are phase transitions at finite T > 0. Since 
the symmetry is now reduced to 0{2) x 0(2) these are, 
generically Kosterlitz-Thouless (KT) transitions. At a 
fixed value of the NLSM coupling go and anisotropy h 
different sequences of phase transitions take place. Some 
details of the phase diagram depends on whether the T — 
quantum phase transition is described by a DFP or a 
BFP. In general there are three situations, shown in the 
phase diagrams of Figl2] (a-c) . 

In Fig. [5] (a) we depict the case in which at T = 
the isotropic 0(4) NLSM is quantum disordered, go > gc, 
corresponding to line c in the T ~ phase diagram of 
Fig 111 For some range of anisotropy h the system re- 
mains disordered at all temperatures. At some critical 
anisotropy /i(5o) there is a T = quantum phase tran- 
sition to a SC or a CDW phase (depending on the sign 
of h), which is in the universality of the 3D XY classical 
model. For T > this transition becomes a KT tran- 
sition and the SC (or the CDW phase depending of the 
case) has power law correlations and not long range or- 
der. Here we assumed that either the T = tetracritical 
point is a DFP or, in the case of of a BFP, that go is 
above the region of the SC-I-CDW coexistence phase. 

In FigO (b) we depict the case in which go = gc (the 
QCP of the 0(4) NLSM), line b in Figd] T^ = only 
for h = 0. The situation depicted in the figure assumes 
that the T = tetracritical point is again of the DFP 
type and hence, that as h increases (in magnitude) one 
does not enter the coexistence phase, but instead the SC 
phase (or CDW depending of the case) which undergoes 
a KT transition to the high temperature phase. 
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Finally, in Figl2] (c) we have the case go < gc, line 
a in Fig HI Once again we have Tc = at /i = as we 
are in 2D. As a function of h (at fixed T) one enters first 
the SC (or CDW) phase (with power law correlations, 
once again). There is a range of h in which the coexis- 
tence phase appears (with KT transitions at both ends) 
and eventually a KT transition to the high temperature 
phase. 

Let us finally discuss, briefly, the finite temperature 
behavior of the 3D stripe state. The main change com- 
pared to the 2D case is that the 0(4) symmetric NLSM 
has now a finite temperature transition, in the universal- 
ity class of the 3D 0(4) Wilson-Fisher fixed point. Hence, 
although Tc for the 0(4) NLSM is lower than that of the 
system with a lower 0{2) x 0(2) symmetry, it is not sup- 
pressed down to T = as in 2D. The 3D case has a phase 
diagram with a more conventionally looking tetracritical 
point. 



V. DISCUSSION 

In this paper we have considered the problem of com- 
petition and coexistence of SC and CDW orders in stripe 
models of strongly correlated systems (e.g. weakly cou- 
pled ladders in the spin gap regime) . We have shown that 
the natural effective field theory in the 2 -I- 1-dimensional 
regime is a spatially anisotropic 0(4) NLSM with ad- 
ditional interactions that break the symmetry down to 
C/(l) X t/(l). We examined the quantum and thermal 
critical behaviors of this system. In the quantum regime 
we used the (/)^ version of this theory (i.e. the Landau- 
Ginzburg- Wilson theory) with 0(4) global symmetry ex- 
plicitly broken down to U{1) x U{1). Using relatively re- 
cent high precision five loop results (improved with Borel- 
Pade resummation of the e expansion) by Calabrese et 
al^ we showed that there is a phase in which SC and 
CDW orders coexist and that this phase transition is 
controlled by a tetra- critical point (likely of the decou- 
pled type). Our results for the quantum and thermal 
phase diagrams (and phase transitions) are summarized 
in the Figures [I] [51 and [3] It is worth noting that our 
RG phase diagrams are essentially the same irrespective 
of the sign of the CDW and SC couplings: we can map 
one problem onto the other by translating the {9i,(j)i) 
fields in equation (j2.2p by tt on all odd-numbered stripes. 
Thus this analysis also describes quantum phase transi- 
tions in a spin-gapped version of the pair-density-wave 
(PDW) state, a phase in which charge and SC orders are 



intertwined,^i25i44^ 

We end by noting that our results are also applicable 
to the problem of a two dimensional array of weakly cou- 
pled spin-1/2 Heisenberg antiferromagnetic spin chains. 
As shown by Senthil and Fisher'^i, a super-spin represen- 
tation used by SF for the 0(4) NLSM describing Neel 
and VBS order parameters for this problem is identical 
to that employed here for us for the problem of coupled 
stripes. The SC and CDW couplings in the stripe prob- 
lem correspond to couplings between the Neel and VBS 
order parameters respectively on nearest neighbor spin 
chains, while departures from Kc = 1 in our problem 
of stripes correspond to departures from the Heisenberg 
point to the XXZ model. Thus, our results suggest that 
in the problem of coupled spin-1/2 chains, there is also 
a phase in which Neel and columnar VBS orders coexist, 
ending at a tetracritical point. Within the coexistence 
region, the VBS singlet order must have spin triplet exci- 
tations. Further, the spatially isotropic 0(4) NLSM the- 
ory of Refs. |39ll40| for the competition of Neel and VBS 
orders in the case of the 2D spin-1/2 Heisenberg antifer- 
romagnetic system admits a deconfined quantum critical 
fixed point at which an anisotropy in the nearest neighbor 
Neel and VBS couplings is a relevant operator; the UV 
stable fixed point reached is the anisotropic 0(4) NLSM 
studied here. It is important to note, however, that our 
0(4) symmetric theory does not have deconfined topolog- 
ical excitations. As we have seen, for {JsCi Jcdw) > 0, 
the competing orders are uniform SC and CDW respec- 
tively, while for (JscJcdw) < 0,the orders are stag- 
gered SC (period-2 pair density wave state) and CDW 
respectively. In the equivalent spin problem, these two 
cases correspond to the competition between antiferro- 
magnetic Neel vs. staggered columnar VBS orders and 
ferromagnetic vs. uniform columnar VBS orders respec- 
tively. 
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